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The function of a real network depends not only on the reliability of its own components, but is 
affected also by the simultaneous operation of other real networks coupled with it. Robustness of 
systems composed of interdependent network layers has been extensively studied in recent years. 
However, the theoretical frameworks developed so far apply only to special models in the limit 
of infinite sizes. These methods are therefore of little help in practical contexts, given that real 
interconnected networks have finite size and their structures are generally not compatible with those 
of graph toy models. Here, we introduce a theoretical method that takes as inputs the adjacency 
matrices of the layers to draw the entire phase diagram for the interconnected network, without the 
need of actually simulating any percolation process. We demonstrate that percolation transitions in 
arbitrary interdependent networks can be understood by decomposing these system into uncoupled 
graphs: the intersection among the layers, and the remainders of the layers. When the intersection 
dominates the remainders, an interconnected network undergoes a continuous percolation transition. 
Conversely, if the intersection is dominated by the contribution of the remainders, the transition 
becomes abrupt even in systems of finite size. We provide examples of real systems that have 
developed interdependent networks sharing a core of “high quality” edges to prevent catastrophic 
failures. 


Percolation is among the most studied topics in statis¬ 
tical physics pQ. The model used to mimic percolation 
processes assumes the existence of an underlying network 
of arbitrary structure. Regular grids are traditionally 
considered to model percolation in materials [21 [3] ■ Com¬ 
plex graphs are instead assumed as underlying supports 
in the analysis of spreading phenomena in social envi¬ 
ronments EQEI, or in robustness studies of technological 
and infrastructural systems EHE]. Once the network has 
been specified, a configuration of the percolation model 
is generated assuming nodes (or sites) present with prob¬ 
ability p. For p = 0, only a disconnected configuration is 
possible. For p = 1 instead, all nodes are within the same 
connected cluster. As the occupation probability varies, 
the network undergoes a structural transition between 
these two extreme configurations. Although there are 
special substrates, e.g., one-dimensional lattices, where 
the percolation transition may be discontinuous, in the 
majority of the cases, random percolation models give 
rise to continuous structural changes [9]. This means 
that the size of the largest cluster in the network, used 
as a proxy for the connectedness of the system, increases 
from the non-percolating to the percolating phases in a 
smooth fashion. 

The percolation transition may become discontinuous 
in a slightly different model involving not just a single 
network, but a system composed of two or more interde¬ 
pendent graphs m- This is a very realistic scenario con¬ 
sidering that many, if not all, real graphs are “coupled” 
with other real networks. Examples can be found in sev¬ 
eral domains: social networks (e.g., Facebook, Twitter, 
etc.) are coupled because they share the same actors HU; 
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Figure 1: Decomposition of interconnected networks into un¬ 
coupled graphs, a) Schematic example of two coupled net¬ 
works A and B. In this representation, nodes of the same 
color are one-to-one interdependent, b) In the percolation 
model, the interconnected system is equivalent to a set of 
three graphs that do not share any edge: the remainders of 
the network layers A and B, and their intersection. 


multimodal transportation networks are composed of dif¬ 
ferent layers (e.g., bus, subway, etc.) that share the same 
locations muni; the functioning of communication and 
power grid systems depend one on the other [TO.. In 
the simplest case, one considers an interconnected sys¬ 
tem composed of only two network layers. Nodes in both 
layers are uniquely identified, but the way they are con¬ 
nected to the other vertices is not necessarily identical 
(see Fig. 1). In the percolation model defined on this 
system, nodes are still present with probability p. Since 
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the networks are interdependent, the presence of a node 
in one layer implies the presence of the same vertex in 
the other layer. However, as p varies, one node may not 
be simultaneously within the largest clusters of both lay¬ 
ers. In such a case, the vertex is said to be outside the 
largest cluster of mutually connected nodes. This is a set 
of nodes identified in a recursive manner, and composed 
of vertices that are simultaneously in the largest clus¬ 
ters of both network layers thanks only to connections 
with other nodes within the set. It has been proved that, 
in infinitely large interconnected systems composed of 
two uncorrelated random networks, the percolation tran¬ 
sition, monitored through the size of the largest cluster of 
mutually connected nodes, is discontinuous [TUI H4l UK] , 
This result has been however shown to not apply to more 
general network models that account for degree correla¬ 
tions ||T6j [T7|. Unfortunately, all these theoretical ap¬ 
proaches have been developed under two special, and 
unrealistic, assumptions. First, they hypothesize that 
network layers are generated according to some kind of 
graph toy model whose topology is not specified by a one- 
zero adjacency matrix, but rather a list of probabilities 
for pairs of nodes to be connected. Second, they apply 
only to the case of infinitely large systems. Real inter¬ 
dependent networks, on the other hand, are composed 
of layers very different from those that can be generated 
with toy models, and they clearly have finite size. In this 
paper, we introduce a novel theoretical approach of di¬ 
rect applicability to the study of percolation transitions 
in real interdependent networks. 

To illustrate our methodology, we start from a simpli¬ 
fied version of a recent theory developed for percolation 
in real isolated networks [181 - 120] . Consider an undirected 
and unweighted graph composed of N nodes and E edges. 
The structure of the graph is encoded by the adjacency 
matrix A 1 a symmetric N x N matrix whose generic ele¬ 
ment Aj :i differs from zero and equals one only if vertices 
i and j share an edge. Without loss of generality, we 
assume that, when all nodes are present in the network, 
the graph is formed by a single connected component. 
Let us consider an arbitrary value of the site occupation 
probability p £ (0,1), and indicate with Si the probabil¬ 
ity that the generic node i is part of the largest cluster. 
The order parameter of the percolation transition is sim¬ 
ply defined as the average of these probabilities over all 
nodes in the graph, i.e., P= A s i■ Note that Si is a 
function of p, but, in the following, we omit this depen¬ 
dence for shortness of notation. As a first attempt, we 
can say that the probability s, for node i to be part of 
the largest cluster is given by 

Si=p[l~ 0 (1 - Sj)] , (1) 

jeW 

where A/] is the set of neighbors of vertex i. The probabil¬ 
ity Si is written as the product of two contributions: (i) 
the probability that the node is occupied; (ii) the prob¬ 
ability that at least one of its neighbors is part of the 
largest cluster. The attempt of Eq. ([2| relies on the so- 


called locally tree-like approximation [9]. In this ansatz, 
neighbors of node i are not directly connected, and this 
allows us to consider the probabilities Sj as independent 
variables. This approximation typically holds in real net¬ 
works m , but does not apply to regular lattices. Intro¬ 
ducing the vectors u and q, whose i-th components are 
respectively rq = ln(l — Si ) and qi = ln(l — Sj/p), we 
can write the set of coupled equations ([I]) into the single 
vectorial equation 

q = Ail (2) 

A trivial solution of Eq. |2]) is given by the configuration 
u = q = 0, corresponding to s = 0 or Si = 0 for all i = 
1,..., N. In the proximity of this configuration, we can 
make use of the truncated Taylor expansion In (1 — x) = 
—x, and Eq. (|2j) reads as 

s = pAs , (3) 

thus an eigenvalue/eigenvector equation. By the Perron- 
Frobenius theorem, the only solution having a physical 
meaning of this equation is obtained by setting p = 1/A 
and s = l, with (A, l) principal eigenpair of the adjacency 
matrix A. This tells us that the solution of Eq. 0 is Si = 
0, for all * = 1,..., ./V, if the site occupation probability 
is smaller than 1/A. In this region, the network is in the 
non-percolating regime. Slightly on the right of 1/A, the 
vector of probabilities s starts to grow in the direction 
of the principal eigenvector of the adjacency matrix, and 
the order parameter is not longer zero. For any value 
of the site occupation probability larger than 1/A, the 
network is in the percolating phase. The site percolation 
threshold obtained using the approximation of Eq. (|T|) is 
thus p c = 1/A. Further, Eq. ([l]) can be solved numerically 
to draw the percolation diagram of the network in this 
approximation. 

The most serious limitation of Eq. (|T|) is to introduce 
a positive feedback among probabilities. An increment 
in the probability Sj produces an increase in the proba¬ 
bilities Sj of the neighbors, which in turn causes an in¬ 
crement in the probability s-i , and so on. To avoid the 
presence of this self-reinforcement effect, we can rewrite 
Eq. (jT]) as 

Si=p[l- 0 (1 ■ (4) 

jeUi 

This equation still relies on the locally tree-like ansatz. 
Here, r^j stands for the probability that node j is part 
of the largest cluster independently of vertex i. We note 
that, while this quantity can be defined for any pair of 
nodes, only contributions given by adjacent vertices play 
a role in Eq. Q. We can think r^j as one of the 2 E 
components of a vector f. In the definition of r, every 
edge (i, j) in the graph is responsible for two entries, 
i.e., Ti^.j and r^-^. For consistency, the probability r^j 
obeys 

r i ^ j =p[l- 0 (I-Tj-hO] • (5) 

fceA/jUd 
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The product of the r.h.s. of the last equation runs 
over all neighbors of node j excluding vertex i. It 
is convenient to rewrite Eq. ([5]) as In (1 — r^j/p) = 
J2 k A jk In (1 - r^k) ~ Aji In (I - rj-n). Defining the 
vectors w and t such that their (i —> j)-th compo¬ 
nents are respectively given by wt^j = ln(l — r^j) and 
Zi^..j = ln(l — ri^.j/p), the system of equations (pi) be¬ 
comes equivalent to the vectorial equation 

t = Mw. (6) 

The generic element of the 2 E x 2 E square matrix M is 
given by 

Mi^j,k—iC ~ ^j,k (f > (f) 

where S X:V is the Kronecker delta function defined as 
fi x ,y = 1 if x = y, and 5 x , y = 0, otherwise. Thus, the 
generic entry of the matrix M is different from zero only 
if the ending node of the edge i —>• j corresponds to the 
starting vertex of the edge k —> t, but the starting and 
ending nodes i and l are different. This matrix is known 
as the non-backtracking matrix of the graph [211 [ 22 ] . 
A trivial solution of the preceding equation is given by 
r = 0, which in turn leads to s = 0. In proximity of 
this configuration, we can still make use of the truncated 
Taylor expansion of the logarithm, and rewrite Eqs. Q 
and © respectively as 

Si = Aijn^j and f=pMf. ( 8 ) 
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Using arguments similar to those applied to Eq. ([3]), we 
can say that, according to Eq. (|s|) , the percolation thresh¬ 
old equals p c = 1/p, with p principal eigenvalue of the 
non-backtracking matrix of the graph, and that slightly 
on the right of the critical point the probability s t grows 
linearly with the sum of the components of the principal 
eigenvector of the non-backtracking matrix correspond¬ 
ing to edges pointing out from node i. The entire per¬ 
colation diagram can be instead obtained by numerically 
solving the system of Eqs. Q and ([5]). 

To summarize, the results presented so far tell us two 
main interesting things. First, the difference between 
the two approaches resides only in the inclusion or exclu¬ 
sion of self-reinforcement effects among local variables. 
In this sense, Eqs. Q and © represent an improve¬ 
ment to Eq. |I]), but both approaches are based on the 
same principles and approximations. This first observa¬ 
tion serves to reunite recent predictions on percolation 
thresholds under the same theory [TM2D] . Second, the 
way in which individual probabilities behave slightly on 
the right of the critical point allow us to understand why 
the prediction of the percolation threshold of Eq. (Jsj) may 
become inaccurate in networks with localized eigenstates 
of the non-backtracking matrix [23]. 

Next, we propose the generalization of the previous 
equations to describe percolation transitions in two in¬ 
terdependent networks. Indicate with A and B the ad¬ 
jacency matrices of the two network layers. Our first 


attempt to write the the probability s, that node i is in 
the largest mutually connected cluster of the system is 
given by 

Si = p [SUB; + (1 - SUbJ SA-Bi Sl3-Ai] > (9) 

where Sx = 1 — YljexO- ~ s i) ' s probability that at 
least one of the nodes j in the set X is part of the largest 
cluster (for the empty set 0, we have S% =0). In the defi¬ 
nition of Eq. ([9]), we have implicitly defined three disjoint 
sets of nodes: ABi = J\fAr\Af B is the set of nodes that are 
neighbors of vertex i in both layers, A — By = \ ABi 

is the set of nodes connected to vertex i only in layer A 
but not in B , and B — Ai = Af B \ ABi is the set of nodes 
that are neighbors of vertex i in layer B but not in A. 
Eq. © essentially states that, given that the vertex is oc¬ 
cupied, the probability s,; for node i of being part of the 
largest mutually connected cluster is given by the sum of 
two contributions: (i) the probability to be connected to 
the largest cluster thanks to at least one vertex that is 
connected to i in both layers; (ii) if the latter condition is 
not true, the probability that node i is connected to the 
largest cluster through at least one node k in layer A and 
one node l in layer B, with k 7 ^ l. Note that, if the net¬ 
work layers are identical, then Eq. ([ 9 ]) correctly reduces 
to Eq. ([I]). In other terms, one can split the set of edges in 
the system in three different subset, and then construct 
three different graphs on the basis of this unique division 
(see Fig. 1): the intersection graph with adjacency ma¬ 
trix given by the Hadamard product of the matrices A 
and B [i.e., the (i,j)~ th element of the adjacency matrix 
is AijBij]; the remnant of network A, where edges be¬ 
tween nodes i and j are present only if A (J (1 — Bij) = 1; 
the remainder of graph B, whose (i,j)-th adjacency ma¬ 
trix element equals Bij(l — A, ? ). If we make use of the 
vector u previously defined, we can write SaBi = 1 — 
ex P E j AijBijUj ], S A -b i = 1 - exp [£N A^( 1 - B^Uj] 
and S^-Ai = 1 — exp [JT Bij( 1 — Aij)uj]. Thus, the nu¬ 
merical solution of Eq. (I9| can be obtained in a certain 
number of iterations, each having a computational com¬ 
plexity that grows at maximum as the number of edges 
present in the denser layer. Unfortunately, the Taylor 
expansion of the r.h.s. of Eq. © gives us only some in¬ 
sights about the structure of the solution, but it does not 
allow to reduce the original problem to a simple eigen¬ 
value/eigenvector equation as in the case of isolated net¬ 
works (see Appendix). 

Also here, we can avoid the presence of self-reinforcing 
mechanisms among variables by excluding already visited 
edges. The equations read as 

Si = p [R-ABi + (1 - B-ABi ) FiA-Bi B-B—Ail , (10) 

and 

=P [l?^B J \{i} + (l _ ^B J \{j}) -^U-BjUi} B B _A j \{i}] , 

(11) 

with Rx t = 1 II j£X (1 U-d) , and the three sets 
ABi, A — Bi and B — Ai are defined as above. If the 


4 




0.0 . . .. . 

0.0 0.5 1.0 

V 



Figure 2: Percolation transition in artificial interconnected networks, a) Percolation diagrams for interdependent Erdos-Renyi 
graphs. We generate network layers with the same average degree (k), and compare results of numerical simulations (large 
symbols) with the solutions of our equations (small symbols). Results are obtained on a single instance of the network model, 
where both layers have size N = 10 4 . Different colors and symbols refer to different values of the average degree (k). b) and 
c) For a given size of the network layers, we generate several instances of the model, and compute the percolation threshold p c 
and the height of the order parameter at criticality P ao (p c )- Both quantities are multiplied by the average degree ( k). Points 
refer to average values obtained over several realizations of the graph model, while error bars stand for standard deviations. We 
use the same symbols and colors as those of panel a. Results are compared with the theoretical expected values p c {k) = 2.4554 
and Poo(Pc)(k) = 1.2564 (dashed black lines), d) We generate a single Erdos-Renyi graph with N = 10 4 and (fc) = 3.0, and 
use it as structure for both layers. We then exchange with probability q the label of every node of layer B with a randomly 
selected vertex, e) Percolation diagrams for single instances of interconnected scale-free networks. Each layer is obtained 
by randomly connecting vertices whose degrees obey the distribution P(k) ~ k \ if k € [5, y/N], and P(k) = 0, otherwise. 
Here, N = 10 4 . f) and g) We compute the average values of p c and Poo(p c ) in several realizations of model composed of 
interdependent scale-free networks with size N. Standard deviations have size comparable with those of the symbol sizes. We 
use the same symbol/color scheme as in panel e. Full lines in panel f stand for best estimates of fits of empirical points with 
the function p c (N) = p c + N~ a . The same type of functions are used in panel g to extrapolate the asymptotic value of -Poo(Pc). 
For any value of the degree exponent, we find that the asymptotic values of both quantities are strictly larger than zero (see 
Appendix), h) We generate a graph with N = 10 4 nodes, and degrees extracted from a power-law distribution with exponent 
7 = 2.5 and support [3, \/iV]. We use this network as structure for both layers. We then exchange with probability q the label 
of every node of layer B with a randomly selected vertex. 


network layers are identical, then Eqs. (10) and (11) 
reduce to Eqs. Q and ^ . If we indicate with P Aa> 
the vector whose components are generated by edges 
present in the intersection graph, w^ AB ' > the vector with 
entries of the type = ln(l — f£^p), and M ^ AB ) 

the non-backtracking matrix obtained from the adja¬ 
cency matrix of intersection between layers, we can write 
= 1 - exp [M {AB) ij}( AB ')]. In a similar spirit, 
we can also write Ra-bPP} = 1 — exp [ M^ A ~ B "> uA A ~ B )} 
and Rs-Aj\{i} = 1 ~ ex P uJ( b_a )], where these 

equations are valid only for edges that belong to either 
layer A or layer B. Obtaining a numerical solution of 
Eqs. (10) and (11) by iteration is thus relatively fast, 


since every iteration has a computational complexity at 
maximum equal to twice the number of edges present in 
the denser network layer. This is a great achievement 
given the high complexity of the algorithm necessary to 
draw the phase diagram for the percolation process in 
interdependent networks by means of direct numerical 
simulations |520- 


Phase diagrams obtained through the numerical solu¬ 
tion of Eqs. (10) and (11) reproduce the results of numer¬ 
ical simulations very accurately. In Fig. 2a for example, 
we consider systems composed of two independent Erdos- 
Renyi network models with different values of the average 
degree (k ), where each network layer is generated by con¬ 
necting pairs of vertices with probability (k) /N. A fun- 
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Figure 3: Percolation transition in interdependent biologi¬ 
cal networks, a) Phase diagram for the multilayer H. sapiens 
protein interaction network [251 [25j . Edges in different layers 
represent diverse type of connections among proteins: direct 
interaction, physical association, and colocalization. When 
analyzing a multiplex with two of these layers, we restrict 
our attention only on the set of nodes present in both layers. 
For each of the three systems formed by two interconnected 
networks that we can generate with this data, we draw the 
percolation diagram by means of numerical simulations (large 
symbols) and numerical solution of our equations (small sym¬ 
bols). b) Phase diagram for the multilayer network of the 
C. elegans connectome [125]. Edges in different layers repre¬ 
sent different types of synaptic junctions among the neurons: 
electrical, chemical monadic, and chemical polyadic. c) De¬ 
composition of the multilayer C. elegans connectome. Rem¬ 
nant of the layer corresponding to electrical junctions, d) in¬ 
tersection among the layers corresponding to electrical and 
chemical monadic interactions and e) remainder of the layer 
corresponding to chemical monadic junctions. In the various 
panels, nodes belonging to the largest connected component 
are visualized with red circles. All other nodes are instead 
represented with blue squares. 


damental feature that the diagrams reveal is the presence 
of a sudden jump in the order parameter P^ at a certain 
threshold p c . We stress that our equations predict the ex¬ 
istence of first-order percolation transitions in networks 
of finite size, and not just in the thermodynamic limit. 
As Figs. 2b and c show, the location of the critical point 
p c , and the height of the jump of the order parameter are 
well described by predictions valid for this type of graph 
models in the limit of infinite size [TO) T¥| US]. We argue 
that a sudden jump in the order parameter is present only 
if the contribution of the remainders dominates the im¬ 
portance of the intersection. This condition is certainly 
verified in interdependent Erdos-Renyi graphs, where the 
intersection is composed of a very small number of edges, 
roughly equal to (k)/ 2 , while the number of edges in 
each remnants is proportional to ( k ) N/ 2. Our intuition 
is fully supported by the results of Fig. 2d. Here, we 
control for the weight of the intersection with respect to 
those of the remnants in a simple fashion [25] , The two 
layers are given by exactly the same network structure. 



Figure 4: Percolation transition in interconnected trans¬ 
portation networks, a) The system is obtained by combining 
Delta and American Airlines routes. We consider only US 
domestic flights operated in January, 2014 m , and construct 
an interconnected network where airports are nodes, and con¬ 
nections on the layers are determined by the existence of at 
least a flight between the two locations. In the percolation di¬ 
agram, large red circles are results of numerical simulations, 
whereas small red circles represent the solutions of our equa¬ 
tions. Blue squares represent susceptibility, a measure of the 
fluctuation across realizations of the percolation model, whose 
peak location is often used as a proxy for the identification of 
the critical threshold p c . b) Same as in a, but for the com¬ 
bination of Delta and United flights, c) Same as in a, but 
for the combination of American Airlines and United flights, 
d, e, and f) Intersection graphs for the systems analyzed 
respectively in panels a, b and c. In the various network 
visualizations, nodes belonging to the largest connected com¬ 
ponent are visualized with red circles. All other nodes are 
instead represented with blue squares. 


Indices of interdependent nodes are however shuffled with 
a given probability q. As q grows, the percolation tran¬ 
sition changes its features: we pass from a continuous 
phase transition for small values of q , through a mix be¬ 
tween a second- and a first-order structural change at 
intermediate values of q , to a discontinuous phase tran¬ 
sition for sufficiently large values of q. The same type 
of considerations hold when network layers are scale-free 
random graphs. The transition is always discontinuous if 
the layers are uncorrelated, so that only remainders are 
present (Fig. 2e). This can be viewed by the existence of 
a finite value of the critical threshold p c (Fig. 2f), and a 
jump of non null height of the order parameter at critical¬ 
ity (Fig. 2g). Still, the nature of the phase transition can 
be tuned from first to second order by simply decreas¬ 
ing the density of the intersection relatively to those of 
the remainders (Fig. 2h). Our argument about the de¬ 
pendence of the nature of the transition on the weight 
of the intersection compared to those of the remnants 
may serve to explain why real interdependent networks 
are not exposed to catastrophic failures El- In Fig. 3, 
we draw the percolation diagrams for two interconnected 
systems of interest in the biological sciences: the H. sapi¬ 
ens protein interaction network [251126] , and the C. el¬ 
egans connectome [26] . Both these interconnected sys- 
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terns undergo continuous percolation transitions. Inter¬ 
estingly, this behavior is not caused by the amount of 
redundancy among layers, but rather the “quality” of 
the edges shared across layers. Connections in the inter¬ 
section graph account, in fact, for less than 10% in five 
out of the six interdependent networks analyzed here. It 
seems therefore that these organisms have developed in¬ 
terconnected networks sharing a core of “high quality” 
edges to prevent catastrophic failures. Whereas the ro¬ 
bustness we observe in biological networks can be viewed 
as the result of a selective evolutionary process, one may 
argue that man-made interdependent systems could have 
been instead not perfectly designed to resist to random 


damages of their components. This is indeed what arise 
from the analysis of the multilayer air transportation net¬ 
work within the US (Fig. 4) [2]5jj30j. The system shows 
fragility, with a sudden jump of the order parameter. On 
the other hand, the height of the jump is not as dramatic 
as observed in random uncorrelated graphs. Major air¬ 
ports all belong to the largest connected component of 
the intersection graph, and their connections constitute 
a set of high quality edges that avoid truly catastrophic 
changes in the connectedness of the entire interdepen¬ 
dent system. From these examples, it seems therefore 
that real interdependent networks may be not so fragile 
as previously believed. 
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7 


7 

Pc 

a 

Poo 

P 

2.3 

0.18 

0.34 

0.01 

0.33 

2.7 

0.25 

0.43 

0.07 

0.39 

3.5 

0.33 

0.60 

0.16 

0.54 

100.0 

0.47 

0.95 

0.32 

1.09 


Table SI: In panel f of Fig. 2, we fitted empirical estimations 
of the percolation threshold p c {N) computed at various sys¬ 
tem sizes N with the function p c (N) = p c + N~ a . In Fig. 2g, 
we perform instead the fit Poo(77) = Px, +N~ l3 for the height 
of the order parameter at criticality. Here, we report the val¬ 
ues of the best estimates of p c , a, Poo and /3 for the different 
values of the degree exponent 7 used in the generation of the 
scale-free networks. 


Appendix 


Since Sa-Bi and Sb~a { calculated at s = 0 are zero, 
the first derivatives of the r.h.s. calculated in s = 0 are 
automatically zero. The Taylor expansion of r.h.s. is 
thus 


2 ’Yhj'Yhk S 3 Sk 
2 Thj Thk s j Sk 


S A-b, S B - A . 
d c do 


„ H + o(sf) = 

s =0 

. H + o(s!) 

s=0 


where the second equality is justified by the fact that 
S A -Bi and S B -At are zero at s = 0. Using the definitions 
of Sa-Bi and S B -Ai we have that 


and 


d q 

— ^A-Bi 

d Sj 


s—0 


Aij{ 1 Bij) 


Taylor expansions 



— By(l Aij) , 

s =0 


An alternative way to arrive to the results of Eq. ^ 
is to use the multidimensional Taylor expansion of the 
r.h.s. of Eq. (JT]) around the trivial solution s = 0 as 


where A and B are the adjacency matrices of the network 
layers. In conclusion, we can approximate Eq. © in 
absence of the intersection term as 


[1 rijeAfi 

= J2k Sk 2 


(1-*,•)] 
^[!-n 


jeMi 


(1-Si)] _ s + o( 8 ?) 


s—0 


— 2 j Aij Sj 


Truncated multidimensional Taylor expansions can be 
used also to reduce Eq. ([5]) to Eq. (|8|. The only differ¬ 
ence here is that the derivatives are taken with respect 
to the variables and the expansion is made around 

the trivial solution r = 0. 


Si ~ 0 Ajj(l Bij)\ [y^ Sj Bjj( 1 Aij)] . 

j j 

With straightforward considerations, we can also insert 
the term accounting for the intersection graph, and write 

Si = p AijBijSj + - [yy Sj Aij(l—Bij)] [ T; Sj -By(l— A^)] . 

3 i 3 


When dealing with Eq. |9]), the Taylor expansion 
should be instead extended to at least the second order. 
Let us first imagine that the intersection graph does not 
contain edges, so that Eq. © reads as 

Si = p S A -Bi Ss-Ai ■ 


This last equation gives us some insights about the struc¬ 
ture of the solution, but it does not allow to reduce the 
original problem to a simple eigenvalue/eigenvector equa¬ 
tion as in the case of isolated networks. Similar consid¬ 
erations can be deduced by taking the Taylor expansion 
of Eq. (111. 


















